report  »  AD-A243  727 


1  10-31-91  Pinal  Technical  15  89  -  14  Auer  Qn 

On  the  Origin  and  Control  of  Large  Coherent 
Structures  in  Turbulent  Shear  Flow 

PE  -  61102P 

PR  -  2307 

SA  -  BS 

G  -  AFOSR  -  89  -  0307 

4.  AufiiUUI  1  1  ■  -  -  .  M  .  . . 

k.  Seifert  and  I.  Wygnanski 

MAAWIAMJi  kaMiU  aiA  million 

Department  of  Fluid  Mechanics  and  Heat  Transfer 

Tel  Aviv  University 

P'^t  Aviv  ISRAEL  AH..;  .• 

aiaoat  Nwam 

AFOSR/HA  1  1 

Building  AlO  • :  i  ^ 

Bolling  AFB  DC  20332-6448 

' L-uC  4  1  199 1  V;  1' 

_ _ 

firpo^ 

ii  iiiiiiuniiy  iiym - - ; - 

iU.  MnMIHHMkAiAliUlV  ITAfUitoT 

Approved  for  public  release;  diaeribution  ia 
unliaiced 

ITie  inhibition  of  two-dimensional  Tollmien-Schlichting  waves  by  active  means  is  well  known. 
Surface  deformation,  heat  input  and  mechanical  means  were  used  to  impose  disturbances  on  the 
boundary  layer  which  have  an  opposite  phase  to  the  waves  existing  in  the  flow.  Thus,  transition 
night  be  delayed  by  suppressing  the  amplitude  of  these  waves.  Natural  disturbances  in  boundary 
ayers  start  as  three  dimensional  wave-packets,  because  they  originate  at  surface  imperfections  or  are 
precipitated  by  temporal  disturbances  in  the  incoming  stream.  The  possibilities  of  controlling  such 
listurbances  are  being  currently  investigated.  In  this  context  the  spatial  interaction  among  three 
I  if\  limensional  wave  trains,  emanating  from  discrete  point-sources,  in  a  boundary  layer  are  discussed, 
i Only  a  local  wave  attenuation  is  feasible  by  activating  two  harmonic,  point-source  disturbances 
I  anywhere  in  the  boundary  layer.  This  was  shown  theoretically  for  a  variety  of  locations,  separation 
j  I  distances  and  phase  delays  between  the  two  sources  and  proven  experimentally  when  the  disturbance 
i  sources  were  displaced  along  the  span.  The  spanwise  phase  gradients  of  the  disturbances,  linked  to 
00  the  streamwise  distance  from  their  point-source  present  a  major  obstacle  to  such  a  simple  attenuation 
iOO  scheme.  These  difficulties  cannot  Im  foreseen  by  considering  3-D  perturbations.  Spatially  distributed 
i  Qji  control  mechanisms  are  therefore  required  for  the  purpose  of  delaying  the  amplification  of 
i  I  concentrated  three  dimensional  disturbances. 
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TRANSITION  TO  TURBULENCE  IN  A  LAMINAR  BOUNDARY  LATER 
AND  ITS  RELATION  TO  LARGE  COHERENT  STRUCTURES  .  , 
IN  TURBULENT  FLOW. 
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Observations  of  natural  transition  is  e}q>erinientally  simple  though  unproductive 
method  of  investigating  the  phenomenon  because  of  the  complexity  of  the  interaction 
among  many  processes  that  occur  almost  simultaneously  during  the  latter  stages  of 
transition.  The  theoretical  analysis  leading  to  the  transition  process  has  its  origin  in 
idealized  two  dimensional  wave- like  perturbations  which,  at  best  are  modulated  along 
the  span.  These  models,  when  backed  by  experiments  (e.g.  Klebanoff  et  al.  1962; 
Kachanov  et  al.  1977;  Saric  &  Reynolds  1980)  predict  the  evolution  of  non  linearities 
and  they  provide  the  notion  that  an  interaction  between  streamwise  and  spanwise  trav¬ 
elling  waves  is  one  of  the  prime  causes  of  transition. 

A  different  experimental  approach  was  tried  out  by  Caster  &  Grant  (1975)  on  the  sup¬ 
position  that  transition  originates  from  a  momentary  disturbance  which  is  hi^ly  con¬ 
centrated  in  a  small  space  namely  at  a  point.  This  supposition  may  stress  the  non 
linear  aspects  of  the  problem  because,  even  a  low-energy  disturbance  originating  by  a 
short  impulse  at  a  point,  might  contain  sufficiently  high  amplitudes  of  various  Fourier 
components  so  as  to  bypass  the  Tollmien-Schlichting  process  and  interact  instantly 
with  the  surrounding  boundary  layer  by  destabilizing  it  and  contaminating  it  with  tur¬ 
bulence.  We  think  that  this  is  the  manner  in  vdiich  turbulent  spots  originate  and  spread 
in  a  laminar  boundary  layer. 


We  have  been  studying  the  spot  over  the  past  IS  years  and  related  its  rate  of  spread 
and  internal  structure  to  known  stability  concepts.  The  relationship  between  the  shape 
of  the  spot  and  its  size  to  the  secondary  breakdown  of  the  Tollmien-Schlichting  wave- 
packet  generated  by  the  spot  was  proven  (Glezer,  Katz  and  Wygnanski  1989).  The 
effect  of  pressure  gradient  on  the  rate  of  spread  of  the  spot  and  on  the  shape  of  its 
boundaries  was  established  (Katz,  Seifert  and  Wygnanski  1990)  and  it  turned  out  that 
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favorable  pressure  gradient  has  a  major  impact  on  the  breakdown  process.  We  are  in  a 
possession  of  detailed  data  pertaining  to  the  internal  structure  of  the  spot  which  has 
still  to  be  properly  analized  (A.  Seifert's  M.Sc.  thesis  was  not  published).  The  persev- 
erence  of  the  the  turbulent  spot  in  a  turbulent  boundary  layer  and  its  probable  connec¬ 
tion  to  the  large  coherent  structures  prevailing  in  the  outer  region  of  a  turbulent  boun¬ 
dary  layer  was  shown  in  the  seventies  (Zilberman  Wygnanski  and  Kaplan  1977)  but 
since  that  time  we  had  been  preoccupied  with  charting  the  structure  of  the  spot  in  a 
laminar  environment.  Consequently  the  interaction  of  the  spot  with  the  turbulence  sur¬ 
rounding  it,  was  not  pursued  in  recent  years  in  spite  of  the  apparent  usefulness  of  this 
approach. 

It  is  quite  possible  that  the  transition  process  is  accelerated  by  the  interaction  of  dis¬ 
turbances  originating  from  different  point-sources  in  the  boundary  layer.  In  order  to 
assess  the  severity  of  these  interactions  and  separate  their  effect  from  other  interac¬ 
tions  resulting  from  either  the  three  dimensionality  of  the  disturbances  or  from  their 
finite  amplitude  the  following  three  part  experiment  was  undertaken.  In  it  only  the 
nature  of  the  interactions  resulting  from  the  superposition  of  two,  small  amplitude  dis¬ 
turbances  were  considered.  These  disturbances  emanated  from  separate  '^int- 
sources"  located  at  different  spanwise  coordinates  on  a  flat  plate.  The  spatial  interac¬ 
tions  between  two  continuous  wave  trains  as  well  as  between  two  wave  packets  were 
investipted  and  special  attention  was  paid  to  the  linear  aspects  of  this  process.  Any 
deviations  from  linearity  were  carefully  scrutinized. 

A  study  of  the  evolution  of  a  disturbance  emanating  from  a  single  harmonic  point- 
source  in  a  Blasius  boundary  layer,  was  initially  repeated  in  order  to  confirm  and 
expand  the  available  data  base  and  provide  a  quick  comparison  for  further  studies. 
These  measurements  were  performed  at  Reynolds  numbers  varying  between  8S0  < 
Rj*  <  1200.  The  three-dimensional,  harmonic  disturbance  is  amenable  to  a  fairly 
simple  linear  analysis  provided  that  the  background  flow  is  parallel  and  the  amplitude 
of  the  disturbance  is  small.  The  eigen  functions  measured  directly  downstream  of  the 
'^int  source”  are  essentially  two-dimensional  and  agree  very  well  with  the  calculated 
eigen-functions  for  this  case.  The  spanwise  distribution  of  the  maximum  amplitude  of 
the  most  prominent  wave  (i.e.  a  wave  whose  dimensionless  frequency  was 
F= 104*10"^  at  990  <R^*  <200)  has  a  saddle  point  directly  downstream  of  the  source 
and  two  maxima  located  along  the  rays  of  Z/X=0.08.  The  double  hump  appearance  of 
the  spanwise  amplitude  distribution  was  shown  to  depend  on  the  frequency  and  on  the 


-  3- 


Reynolds  number.  The  present  results  are  in  substantial  agreement  with  the  measure¬ 
ments  of  Kachanov  and  with  the  linear  calculations  of  Mack.  The  shape  of  the  mea¬ 
sured  amplitudes  and  phases  of  the  streamwise  velocity  perturbations  changes  consid¬ 
erably  along  the  span  due  to  the  coupling  between  the  Orr-Sommerfeld  modes  and 
vertical  vorticity  modes.  The  measurements  are  in  good  agreement  with  the  calculated 
linear,  three-dimensional  eigen  functions.  The  harmonic  point  source  disturbance 
evolves  Li  space  according  to  the  linear  model  provided  the  forcing  amplitude  does  not 
exceed  3%  within  the  range  of  R^*  considered,  or  less  than  1%  far  downstream  from 
the  source. 

The  interactions  of  two  harmonic  wave-trains  emanating  with  the  same  phase,  ampli¬ 
tude  and  frequency  from  two  sources  which  are  only  displaced  along  the  span  were  in- 
vestipted  in  the  second  stage  of  the  experiment.  Surveys  of  velocity  fluctuations  were 
made  at  various  distances  from  both  sources.  The  data  for  each  survey  was  recorded 
three  times:  once  while  both  harmonic  disturbances  were  simultaneously  activated  and 
twice  when  only  one  of  the  disturbances  was  alternately  activated.  Only  in  this 
manner  one  could  asses  the  validity  of  the  linear  superposition  of  these  perturbations. 
At  sufficiently  low  amplitudes  of  the  initial  harmonic  perturbation,  a  complex  summa¬ 
tion  of  the  waves  emanating  from  each  source  separately,  represented  quite  well  the 
flow  pattern  observed  when  both  sources  operated  concurrently. 

Wave  cancellation  by  two  harmonic  point-source  disturbances  originating  from  dif¬ 
ferent  locations  on  the  surface  (i.e.  displaced  in  both  x  and  z)  are  virtually  impossible 
due  to  the  steep  spanwise  phase  gradients  of  the  waves  further  downstream.  This  type 
of  difficulty  was  not  encountered  whenever  the  perturbations  were  two-dimensional. 
Consequently,  much  more  sophisticated  control  mechanisms  are  required  in  order  to 
delay  the  amplification  of  three  dimensional  disturbances  in  the  direction  of  streaming 
than  hitherto  contemplated. 

In  the  second  stage  ''f  this  experiment,  the  continuous  oscillatory  perturbation  was 
replaced  by  a  pulse,  which  in  turn,  generated  a  wave  packet  emanating  from  a  point 
source.  The  excitation  amplitudes  were  low  enough  so  that  the  streamwise  evolution 
of  the  disturbance  was  essentially  linear  over  a  considerable  distance.  The  amplitude 
and  phase  of  each  Fourier  component  of  the  perturbation  in  the  packet  was  identical 
to  the  corresponding  perturbation  in  the  wave  train,  even  the  appearance  of  the 
double- hump  in  the  spanwise  distribution  of  amplitudes  reoccurred  in  the  wave  packet 
for  the  predominant  frequency  observed.  The  flow  field  associated  with  the  packet 
was  not  only  mapped  outside  the  boundary  layer,  as  it  was  done  by  Caster,  but  all 
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across  the  boundary  layer  and  in  particular  near  the  critical  layer  (near  the  surface) 
where  the  amplitudes  are  higji. 

When  the  spanwise  interaction  of  two  wave  packets  was  considered,  a  weakly  non 
linear  interaction  was  noticed  at  large  values  of  R^*  for  the  lowest  wave  numbers 
considered.  The  coexistence  of  a  broad  spectrum  of  velocity  fluctuations  in  the  inter¬ 
acting  wave  packets  may  produce  a  large  number  of  secondary  interactions  which  are 
not  feasible  during  the  interaction  of  two  harmonic  wave  trains.  It  appears  that  the 
subharmonic  form  of  interaction  emerges  as  being  the  leading  one  in  the  range  of  par¬ 
ameters  considered. 
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Abstract 


The  inhibition  of  two-dimensional  Tollniien- 
Schiichting  waves  by  active  means  is  well 
known.  Surface  deformation,  heat  input  and 
mechanical  means  were  used  in'  impose  distur¬ 
bances  on  the  boundary  layer  which  have  an 
opposite  pliase  to  the  waves  existing  in  the  flow. 
Thus,  provided  that  two-dimensional  T-S  waves 
lead  to  turbulence,  transition  might  be  delayed 
by  suppressing  the  amplitude  of  these  waves. 
Natural  disturbances  in  boundary  layers  are 
manifested  as  tliree  dimensional  wave-packets, 
because  tliey  originate  at  surface  imperfections 
and  are  precipitated  by  temporal  disturbances  in 
the  incoming  stream.  The  possibilities  of  con¬ 
trolling  such  disturbances  are  being  currently 
investigated.  In  ’.his  context  the  spatial  interac¬ 
tion  among  three  dimensional  wave  trains,  ema¬ 
nating  from  discrete  point-sources,  in  a  Blasius 
boundary  layer  are  discussed. 

The  evolution  of  an  isolated  harmonic  distur¬ 
bance  emanating  from  a  point-source  was  in¬ 
itially  considered  in  order  to  provide  a  reference 
for  tlie  mure  complex  configurations.  The 
experiments  were  done  at  Reynolds  numbers 
varying  between  820  <  R^«  <  1150  and  the  data 

was  compared  to  linear  calculations  assuming 
that  the  background  flow  is  parallel  and  the 
amplitude  of  the  disturbance  is  small.  The  mea¬ 
surements  validated  tite  linear  calculations 
within  the  range  of  amplitudes  and  R^«  consi¬ 
dered.  They  are  also  in  substantial  agreement 
with  the  measurements  of  Kachanov  (Ref.  1) 
and  with  the  predictions  of  Mack  (Ref.  2). 


The  investigation  of  the  interaction  between  two 
harmonic  wave-trains  emanating,  with  the  same 
amplitude  and  frequency,  from  two  point- 
sources  displaced  along  the  span  followed.  The 
validity  of  the  linear  superposition  of  these  dis¬ 
turbances  was  confirmed  experimentally 
because,  at  sufficiently  low  amplitudes  of  the 
initial  harmonic  perturbation,  a  complex  sum¬ 
mation  of  the  waves  emanating  from  each 
source  separately,  also  represented  the  resulting 
flow  pattern  when  both  sources  operated  con¬ 
currently. 

Only  a  local  wave  attenuation  is  feasible  by  acti¬ 
vating  two  harmonic,  point-source  disturbances 
anywhere  in  the  boundary  layer.  This  was 
shown  theoretically  for  a  variety  of  locations, 
separation  distances  and  phase  delays  between 
the  two  sources  and  proven  experimentally  when 
the  disturbance  sources  were  displaced  along  the 
span.  The  spanwise  phase  gradients  of  the  dis¬ 
turbances,  linked  to  the  streamwise  distance 
from  their  point-source  present  a  major  obstacle 
to  such  a  simple  attenuation  scheme.  These  dif¬ 
ficulties  could  not  be  foreseen  by  considering 
two-dimensional  perturbations.  Spatially  distri¬ 
buted  control  mechanisms  are  therefore  required 
for  the  purpose  of  delaying  the  amplification  of 
concentrated  three  dimensional  disturbances, 

either  steady  or  pulsating,  with  increasing  R^«. 
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List  of  Symbols 


Introduction 


A 

Amplitude,  (%  of  Uqo] 

C 

Phase  velocity,  c=w/a 

f 

frequency,  (HzJ 

F 

dimensionless  frequency,  xlO® 

Uoo 

Reynolds  number,  R^*=(Uoo  6*)/u 

Re^j 

Reynolds  number,  Rejj=(Uoo  X)/t/ 

Re^ 

Reynolds  number,  Rej^=(U<3o  Z)/u 

U 

mean  velocity  in  X  direction 

u 

fluctuating  X  component  of  velocity 

V,v 

Y  component  of  velocity 

W,w 

Z  component  of  velocity 

X 

streamwise  coordinate 

X, 

streamwise  distance  from  disturbance 
source 

Xp 

X  location  of  disturbance  source 

Y 

coordinate  normal  to  the  wall 

^nia 

Y  of  disturbance  maximum  amplitude 

Z 

spanwise  coordinate 

Zs 

spanwise  distance  from  disturbance 
source,  Z3=Z-Zp 

Zp 

Z  location  of  disturbance  source 

Greek 

• 

‘»r 

real  part  of  streamwise  wavenumber, 

af=2x/Aj( 

Qi 

amplification  factor 

spanwise  wavenumber,  real, 

6* 

boundary  layer  displacement  thickness 

h 

vertical  component  of  vorticity 

K 

3D  wave  number,  <c*=q*+^* 

streamwise  wavelength 

spanwise  wavelength 

V 

kinematic  viscosity 

<A 

phase 

'ii 

stream  functions 

X 

complex  phase 

ill 

circular  frequency,  iir*2ir( 

Abbreviations 

FSP 

Falkner-Skan  Parameter,  usually  labeled 
/3,  /J»2m/(m+U  .  where  U<x,  «  X”' 

HPS 

A  harmonic  disturbance  emanating  from 
an  isolated  Harmonic- Point-Source 

MA 

Maximum  Amplitude 

TUPS 

A  harmonic  disturbance  emanating  from 
Two-llarmonic-Point-Sources 

TS 

Tollmein-Schlichting  waves 

2D 

Two  Dimensional 

3D 

Three  Dimensional 

Background 

Velocity  perturbations  occurring  naturally  in  a 
boundary  layer  are  seldom  harmonic  or  two 
dimensional.  They  often  originate  at  surface 
imperfections  and  are  precipitated  by  temporally 
random  disturbances  in  the  free  stream.  Nev¬ 
ertheless,  the  most  detailed  information  available 
about  the  growth  or  decay  of  disturbances  in 
boundary  layers  focuses  on  a  somewhat  artificial 
case  of  a  harmonic  motion  in  two  dimensions. 
The  current  understanding  of  the  possible  means 
of  controlling  such  disturbances  is  even  more 
restricted.  Tollmien-Schlichting  waves  of  identi¬ 
cal  frequency,  emanating  from  two,  2D  har¬ 
monic  sources  displaced  in  the  direction  of  stre¬ 
aming,  may  either  reinforce  or  weaken  one 
another  depending  on  the  phase  relationship 
between  them.  The  possibility  of  wave  cancella¬ 
tion  is  attractive  because  it  may  be  used  to  delay 
the  transition  from  laminar  to  turbulent  flow. 
Two  dimensional  wave  interactions  were  first 
investigated  by  Wehrmann  (Ref.  3),  by  Liep- 
mann  et.  al  (Ref.  4),  Thomas  (Ref.  5)  and  most 
recently  by  Popator  and  Saric  (Ref.  6).  A  single 
mode  of  oblique  waves  could  be  controlled  in  a 
similar  manner. 

A  suppression  of  complex  3D  disturbances  ori¬ 
ginating  randomly  by  a  variety  of  sources  repre¬ 
sents  a  difficult  task.  Isolated  disturbances  ema¬ 
nating  from  a  single  point  source  were  studied 
by  Caster  (Refrences  7  and  8)  and  Muck  (Ref. 
2),  following  some  pioneering  work  of  Benjamin 
(Ref.  9)  and  of  Criminate  and  ICovasznay  (Ref. 

10) .  The  disturbances  investigated  are  either 
continuous  or  time  dependent  and  consequently 
they  were  generated  by  either  a  harmonic  point- 
source  or  by  a  pulse.  Caster  (Ref.  7)  showed 
that  three  dimensional  wave-packets,  constitut¬ 
ing  of  a  finite  band  of  frequencies  and  spanwise 
wave  numbers  propagating  in  a  boundary  layer, 
can  be  generated  by  a  localized  impulsive  dis¬ 
turbance  (e.g.  a  brief  and  a  tiny  vertical  jet). 
Such  a  disturbance  produces  a  flat,  broad-band 
spectrum  which  is  selectively  amplified  by  the 
boundary  layer  to  contain  a  dominant  band  of 
2D  and  3D  wave  numbers.  Caster's  theoretical 
model  was  based  on  the  integration  of  eigen 
modes  calculated  from  the  On -Sommerfeld 
equation  making  use  of  Squire's  theorem  (Ref. 

11) .  These  calculations  resembled  the  observed, 
initial  development  of  the  wave-packet  at  the 
outer  edge  of  the  boundary  layer  (Ref.  8). 
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Mack  and  Kendall  (Ref.  2)  and  Kachanov  (Ref. 
I)  investigated  the  evolution  of  a  wave  train 
produced  by  an  harmonic-point-source  (HPS)  in 
a  Blasius  boundary  layer.  Mack's  calculations  in¬ 
dicated  that  the  maximum  amplitude  of  the  dis¬ 
turbance  lies  on  the  ray  emanating  from  the 
source  and  defined  by  Zj  /  Xj  =  0.05  to  0.09 
(for  F=00  at  R^*=800  to  1000),  rather  than  on 

the  plane  of  symmetry  as  expected  from  the 
simplified  interpretation  of  Squire's  theorem. 
These  calculations  are  in  good  agreement  with 
experimental  observations  (References  I  and  2). 

The  interaction  between  two  wave  trains  ema¬ 
nating  from  spatially  separated,  two  harmonic- 
point-sourccs  located  in  a  fully  developed  Bla¬ 
sius  boundary  layer  is  considered  presently.  It  is 
a  much  simpler  type  of  interaction  than  the  one 
occurring  between  wave  packets  yet  it  contains 
some  important  three  dimensional  ingredients 
which  need  to  be  considered  in  any  attempt  of 
active  control  of  the  transition  process.  This 
simpler  kind  of  interaction  enables  us  to  study 
separately  the  effect  of  three  dimensionality, 
finite  amplitude  and  wide  frequency  content 
which  are  alt  present  simultaneously  in  the  wave 
packet  interaction.  Special  attention  was  paid  to 
the  complex  30  linear  interactions  which  might 
precipitate  large  amplitudes  even  when  the 
amplitudes  of  the  individual,  non  interacting 
perturbations  are  small. 

The  project  was  subdivided  into  three  interme¬ 
diate  steps  consisting  of;  (i)  an  experimental  and 
theoretical  investigation  of  a  single  harmonic 
disturbance  originating  from  a  point;  (ii)  an  in¬ 
teraction  of  two  harmonic  disturbances  dis¬ 
placed  in  the  span  wise  direction  and  (iii)  a 
calculation  simulating  an  interaction  of  two 
harmonic  disturbances  displaced  in  the  stream- 
wise  direction. 


Equation  (2.1)  represents  the  3D  equivalent  of 
the  Orr-Sommerfeld  (OS)  equation  from  which 
both  V  and  the  eigen  values  can  be  calcu¬ 
lated.  The  vertical  vorticity  component  fy  which 
appears  in  equations  (2. 1  -4)  is  given  by: 

'y  dz  dx 

and  therfore  equation  (2.2)  is  referred  to  as  the 
vertical  vorticity  equation.  In  the  spatial  case  a, 
is  complex  while  u  and  are  real. 

The  boundary  conditions  for  the  OS  equation 
are: 

flu 

at  Y«o,oo  (4.1) 

and  for  the  vertical  vorticity  equation: 

fy  »  0  at  Y  -  0  ,  oo  (4.2) 

In  this  context  it  is  worth  noting  that  the  well 
known  Squire  transformation  (Ref.  II)  is  com¬ 
monly  interpreted  to  imply  that  every  3D 
normal  mode  has  an  equivalent  2D  normal 
mode,  existing  at  a  lower  Reynolds  number. 
However,  this  simplified  assumption  is  not 
always  helpful,  since  certain  combinations  of 
frequency,  Reynolds  number  and  wave  orienta¬ 
tion,  may  actually  lead  to  higher  3D  amplifica¬ 
tion  rates. 


Linear  Theory 


Consider  the  3D  disturbance  to  be  of  the  fami¬ 


liar  form: 

“i  (X.Y.Z.t)  =  “j  (Y)  e  '  <®^+^2:-wt)  (I) 


where  j  =  1, 2, 3, 4  corresponds  to  the  three  com¬ 
ponents  of  velocity  u,  v,  w,  and  pressure  fluctu¬ 
ations  respectively  and  substitute  it  into  the  lin¬ 
earized  equations  of  motion.  The  result  is  a  set 
of  equations  derived  in  Ref.  12: 


The  3D  Orr-Somerfeld  equation  was  solved 
numerically  by  using  a  shooting  procedure 
described  in  detail  elsewhere  (Ref.  13,  for  a 
solution  of  the  2D  case).  Proper  changes  have 
been  implemented  in  the  above  procedure  for 
the  integration  of  the  full  3D  equation  (i.e  Eq. 
2.1-4).  For  a  convenient  comparison  with  exper¬ 
iment,  the  complex  eigen  functions  are  repre¬ 
sented  by  an  amplitude  and  a  phase: 
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Numerical  Prediction  of  t!>e  UPS  Maximum 
Amplitude 


The  normal  mode  representation  of  a  three 
dimensional,  small  amplitude,  harmonic  distur¬ 
bance  has  the  following  form  (Ref.  7); 


X-®  “  Exp  . 

i 

I  QdX  +  ^Z  -  wt 

Jxs 

> 

- 

(6) 


where  is  a  factor  that  accounts  for  the 

effect  of  the  boundary  layer  divergence.  Mack 
(Rfi  proposed  a. similar  expression  for  the 
amplitude  of  the  streamwise  component  of  the 
velocity  fluctuations; 


u  (X,Z,t,ai)  =  Exp(iwt) 


Expi^/3,X,Z)d^  (7) 


and  for  the  complex  phase: 

0(X.Z)=  a  (X;^,u)  dX  +  (8) 

JXs 

where  a  is  the  complex  eigen  value  obtained 
from  the  30  OS  equation  and  the  spanwise  loca¬ 
tion  of  the  source  is  at  Zp*0.  This  model,  like 
the  previous  one  (Ref.  7),  also  ignores  the  shape 
of  the  eigen  functions. 


The  numerical  scheme  used  presently  for  calcu¬ 
lating  the  evolution  of  the  iiPS  disturbance  con¬ 
tained  (he  following  steps: 

1)  The  eigen  values  of  (he  spatial  30  OS 
equation  (Eq.  3.4)  were  calculated  for  the 
Olasius  velocity  profile  in  the  range  of  700 
<  R  *  <  1200  at  increments  of  R,*  «  25. 

o  6 

The  inclination  of  (he  wave  crests  relative 
to  the  Z  axis  varied  from  0-70®  at  I®  incre¬ 
ment.  Advantage  was  taken  of  the  anti¬ 
symmetry  in  /}. 

2)  Calculation  of  the  amplitudes  and  the 
phases  over  the  entire  region  considered  in 

the  e.xperinient  i.e.  for  700  <  R^*  <  1200 

and  I  Zj  /  Xj  I  <  0.3. 

Experiment 

The  evolution  of  a  single  wave  train  and  the  in¬ 
teraction  between  two  wave  trains  emanating 
from  separate  sources,  located  in  a  Blasius 
boundary  layer  and  displaced  in  the  spanwise 
direction,  were  investigated  experimentally.  The 
experinient  was  conducted  in  the  close  loop,  low 
turbulence  wind  tunnel  at  Tel-Aviv  university. 


A  fully  automated,  computer  controlled  experi¬ 
ment  was  constructed,  which  encompassed  auto¬ 
matic  hot-wire  calibration  procedure  and  pruL. 
traversing  in  the  proximity  of  the  wall.  A  pol¬ 
ished  aluminum  plate  19  mm  thick  and  2.5  m 
long  was  installed  vertically  in  the  10'  long  rem¬ 
ovable  test  section.  The  plate  has  an  elliptic 
leading  edge  and  a  trailing-edge  flap  in  order  to 
control  the  circulation  around  the  plate  and 
through  it  the  location  of  the  leading-edge  stag¬ 
nation  line.  The  turbulent  corner  flow  originat¬ 
ing  at  the  juncture  of  the  plate  and  the  side 
walls  was  forced  to  bleed  to  the  backward  side 
of  the  plate  through  adjustable  slots.  This  was 
achieved  by  a  slight  increase  in  the  static  pres¬ 
sure  on  the  working  side  of  the  plate  which 
helped  to  maintain  a  long  fetch  of  laminar  flow 
on  that  surface.  The  free  stream  velocity  was 
maintained  at  7.5  m/s  throughout.  The  calcu¬ 
lated  Reynolds  number  based  on  the  measured 
displacement  thickness  agreed  with  the  expected 
values  for  the  Olasius  boundary  layer  as  shown 
in  figure  I. 


Fi»  I  R,*  vs  ifx 


The  spanwise  distribution  of  R^*  in  the  region 

of  interest  was  maintained  at  its  nominal  v^lue 
by  5%.  The  top  and  bottom  walls  of  the  test 
section  were  adjusted  until  the  measured 
Falkner-Skan  Parameter  (FSP)  vanished,  or  at 
least  did  not  exceed  the  value  t0.05  which  was 
the  possible  resolution  of  the  measurement.  The 
FSP  is  a  very  important  parameter  strongly 
affecting  the  critical  Reynolds  number  and  the 
linear  amplification  factor. 

Controlled  harmonic  disturbances  were  gener¬ 
ated  by  two  miniature  hearing-aid  earphones 
embedded  on  the  rear  side  of  the  flat  plate  and 
ejecting  tiny  jets  of  air  through  0.5  mm  holes. 
These  sources  were  'ocated  50  cm  from  the 

leading  edge  of  the  plate,  corresponding  to  R  ♦ 

0 

=  700.  The  spanwise  separation  between  the 
two  sources  was  80  mm  which  was  equivalent 
to  approximately  50  local  displacement 
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i'ig.  2  Flai  plate  apparatus,  spreading  rate  ot  the  two  HPS  disturbances  and  overlap  regions. 


thicknesses  at  the  free  stream  velocity  of  7.5 
m/sec  and  amounted  two  wave  lengths  of  a 
typical  perturbation  at  the  dimensionless  fre¬ 
quency  of  F=>I04.  Doth  sources  were  activated 
by  the  same  input  signal  and  their  output  amp¬ 
litudes  were  matched  in  order  to  produce  equal 
atnplitudes  at  any  desired  phase  shift. 

The  disturbances  from  a  single  harmonic-point- 
source  spread  in  the  spanwise  direction  at  half 
an  included  angle  of  8-10**,  consequently  the 
prescribed  spanwise  separation  between  the 
sources  allowed  each  wave  train  to  develop  in¬ 
dependently  before  it  overlapped  and  interacted 
with  its  neighbor  (see  Figure  2).  The  interaction 
experiment  was  done  in  three  steps.  In  the  first 
step  both  sources  operated  simultaneously  while 
subsequently  each  source  operated  alone  under 
identical  flow  conditions.  In  this  way,  any  resi¬ 
dual  asymmetry  in  the  flow  or  any  uncontroll- 
aole  perturbation  in  the  tunnel  would  affect  all 
three  stages  of  the  experimen’  equally  and 
could  be  subtracted  out.  The  amplitude  and 
phase  distributions  of  the  streamwise  compo¬ 
nent  of  the  velocity  fluctuation  were  calculated 
from  data  sampled  across  the  entire  boundary 
layer  at  every  X,Z  location.  More  than  30 
points  of  the  streamwise  velocity  component 
were  measured  within  the  boundary  layer 
between  the  wall  and  and  the  Y  location  cor¬ 
responding  to  46  .  The  maximum  amplitudes 
presented  were  determined  with  respect  to  Y  at 
every  X,Z  location  while  the  phase  reference 


was  chosen  at  the  location  at  which  the  ampli¬ 
tude  was  maximum  (i.e  Y^g).  The  flow  was 
interrogated  along  the  center-plane  of  every 
isolated  harmonic  disturbance  (i.e  Z^^Zp),  along 
the  center-plane  of  the  combined  disturbance 
(i.e.  at  Z=0)  and  across  the  entire  span  of  the 
disturbance  at  a  few  select  stations.  The 

effect  of  constant  phase  mismaich  between  the 
two  harmonic  disturbances  was  documented  in 
order  to  study  the  outcome  of  such  a  mismatch 
hoping  to  achieve  a  favorable  interaction  and 
possibly  a  suppression  of  this  type  of  distur¬ 
bances.  The  two  harmonic  sources  were  oper¬ 
ated  at  almost  equal  amplitudes  but  at  a  phase 
differences  of  0,  x/2  and  jt. 
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Results  and  Discussion 


The  evolution  of  an  isolated  harmonic  point- 
source  disturbance 

The  shape  of  the  amplitude  and  phase  distribu¬ 
tions  vs  Y  of  the  velocity  fluctuations  change 
along  the  span  of  the  HPS,  implying  that  there 
is  a  concomitant  change  in  the  form  of  the  eigen 
function  along  the  span,  in  contradistinction  to 
the  20  case.  The  amplitude  and  phase  distribu¬ 
tions  measured  directly  downstream  of  the  single 
source  (i.e.  at  Z=Zp)  agreed  very  well  with  the 
eigen  functions  calculated  for  the  purely  20 
case  (Figure  3a j. 

This  is  nut  surprising,  because  symmetry 
requirements  dictate  that  the  dominant  wave  on 
the  center- plane  must  be  two  dimensional  (i.e. 
/3=0),  as  demonstrated  on  the  top  portion  of 
figure  4.  The  spanwise  distribution  of  the  maxi¬ 
mum  am|)litude  and  the  phase  shift  measured  at 
the  Y  location  at  which  the  maximum  amplitude 
occurs  relative  to  its  occurrence  on  the  plane  of 
symmetry  are  plotted  on  the  lower  part  of 
figure  4. 1  he  solid  lines,  in  this  figure,  represent 
the  best  fit  to  the  data.  The  spanwise  wave 
number  becontes  infinitely  large  toward  Zp^O 
and  is  anti-symmetrical  with  respect  to  the 
centerline,  on  which  it  is  suppose  to  vanish 
implying  that  the  disturbance  is  two  dimensional 
there.  Away  from  the  center-plane  (e.g.  on  the 
ray  Zj  /  X5=0.136),  where  the  oblique  waves 
are  dominant  and  /36  =0.2,  the  experimentally 
determined  eigen  functions  are  quite  different 
(Figure  3b).  The  calculated  phase  and  amplitude 
distributions  of  the  u  component  of  the  velocity 
perturbation  for  P6*=Q.2,  F=I04,  and  R^osllOO 

are  in  excellent  agreement  with  the  measured 
data.  Similar  agreement  was  found  at  other  fre¬ 
quencies  and  spatial  locations,  as  it  is  demon¬ 
strated  by  the  results  presented  in  figure  5.  It 
should  be  mentioned  that  the  data  is  normalized 
whereupon  the  maximum  amplitude  and  one 
point  in  the  phase  distribution  are  matched  with 
the  calculated  values.  An  increase  in  the  forcing 
amplitude  did  not  cause  measurable  distortions 
of  the  eigen  functions  or  a  noticeable  change  in 
the  amplification  rate  provided  the  forcing 
amplitude  was  less  than  3%  of  the  free  stream 
velocity  (or  a  local  amplitude  of  1%  further 
downstream). 


A/Ainoi  (ji  [d«gj 


Fig.  3  Measured  and  calculated  eigen  functions 
on  (a)  and  off-centerline  (b),  F=104  at 
R^**990  (Xs=35  cm) 


Fig.  4  Experimental  HPS  spanwise  waven¬ 
umber,  wave  length,  phase  and  maxi¬ 
mum  amplitude  vs  Z,  F>I04  at  R^**990 

(Xs-35cm). 
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rig.  5  Comparison  between  measured  and  caicuiated  eigen  functions  on  and  off-centerline, 
l•=138  at  R,**990. 


Fig.  6  Normalized  maximum  amplitudes  vs 
for  F=69,  104  and  138,  comparison 

6 

between  experiment  and  theory. 


F*I04 


Fig.  7  Calculated  HPS  maximum  amplitudes  vs 
on  Zs=0,  F*I04,  for  Falkner-Skan 

parameters  of  -0.03,  0  and  0.03. 


A  comparison  between  the  measured  and  the 
calculated,  maximum  amplitude  distributions  on 
the  plane  of  symmetry,  for  a  variety  of  dimen¬ 
sionless  frequencies,  is  made  in  figure  6.  The 
calculations  are  based  on  Mack's  model  and  are 
in  fair  agreement  with  the  experimental  obser¬ 
vations. 
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Fig.  8  Comparison  between  experimental  and  theoretical  maximum  amplitude  and  pliase  at 
Y=Yn,a  vs  Z  at  R^*=850  and  990,  F=104  and  138. 


When  assessing  the  quality  of  this  cmnparison 
one  should  consider  the  tremendous  sensitivity 
of  the  measured  amplification  rate  to  pressure 
gradients.  The  stream  wise  amplification  of  a 
dominant  harmonic  disturbance  is  presented  in 
figure  7  for  the  range  of  distances  used  in  the 
present  experiment.  The  three  curves  shown 
correspond  to  different  Falkner-Skan  parame¬ 
ters,  ranging  between  ±0.05  values.  U  should  be 
remembered  that  the  difference  between  a  velo¬ 
city  profile  corresponding  to  a  FSP  of  0.05  and 
of  0  is  hardly  discernible  experimentally,  yet  the 
differences  in  the  integrated  amplification  rate 
over  the  distances  considered  presently  may 
result  in  amplitudes  which  are  orders  of  magni¬ 
tude  apart.  Consequently  the  disparity  between 
theory  and  experiment  shown  in  figure  6  is  in¬ 
significant.  The  sensitivity  of  the  streamwise 
amplification  to  pressure  gradient  dissuaded  us 
from  refining  the  theoretical  model  by  including 
in  it  the  slight,  streamwise  divergence  of  the 
laminar  boundary  layer. 

Spanwise  amplitude  and  phase  distributions, 
corresponding  to  the  Y  location  at  which  the 
amplitude  is  maximum  were  also  calculated 
using  Mack's  model.  A  sample  of  these  calcula¬ 


tions  is  compared  with  experimental  data 
(Figure  8)  for  two  streamwise  locations 
(K,.«»850  and  990)  and  for  two  dimensionless 

frequencies  (Fsl04  and  138). 

It  should  be  noted  that  the  spanwise  amplitude 
and  phase  distribution  depends  on  the  distance 
from  the  source,  on  the  frequency  and  on 
the  streamwise  location  of  the  source  itself.  It  is 
anticipated  therefore,  that  only  a  local  and  par¬ 
tial  suppression  of  the  wave  can  be  achieved  by 
the  destructive  spanwise  interaction  between 
two,  locally  concentrated,  harmonic  distur¬ 
bances. 
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(Xs-55cm). 


Spanwise  interaction 

Three  spanwise  distributions  of  the  maximum 
amplitudes  observed  at  a  given  streamwise  loca¬ 
tion  are  plotted  in  figure  9  together  with  their 
respective  distribution  of  phase-angles.  These 
results  were  obtained  while  the  two  spanwise 
separated  point-sources  of  harmonic  distur¬ 
bances  were  operating  simultaneously,  at  almost 
equal  amplitudes.The  phase  difference  between 
the  sources  was  controlled  and  the  data  pre¬ 
sented  corresponds  to  three  principle  phase 
differences  of  0,  »/2  and  »  radians. 


The  perturbation  resulting  from  such  two  inter¬ 
acting  disturbances  may  be  effectively  recon¬ 
structed  from  the  data  collected  for  each  isolated 
disturbance,  by  a  linear  superposition,  provided 
that  the  shape  of  the  eigen  function  at  every 
X,Z  location  is  known  and  the  spanwise  separa¬ 
tion  and  phase  shift  are  taken  into  consideration. 
In  order  to  form  a  correct  linear  superposition 
of  the  two  disturbances,  one  should  sum  up  the 
two  complex  eigen  functions  acquired  at  the 
same  X,Z  location,  in  the  form 


One  may  observe  from  this  illustration  that  the 
initial  phase  difference  does  not  affect  the  max¬ 
imum  amplitude  attained  at  some  arbitrary  X,Z 
location  in  the  boundary  layer.  The  control  of 
the  phase-shift  between  the  harmonic  distur¬ 
bances  can  only  displace  the  location  at  which 
the  maximum  amplitude  occurs  along  the  span 
but  not  suppress  this  amplitude.  This  is  not 
surprising  when  one  takes  into  consideration  the 
steep  spanwise  phase  gradients  occurring  down¬ 
stream  of  an  isolated,  harmonic-point-source 
disturbance. 


♦sup  (y)  -  A^e*^*  +  A,e‘^*  (9) 

Aj,  Aj,  and  are  functions  of  Y.  The  sup¬ 
erimposed  signal  should  be  compared  with  the 
ensuing  disturbance  generated  by  the  two  har¬ 
monic  point-sources  (TUPS)  operating  in 
unison.  Any  deviation  of  the  measured  ampli¬ 
tude  or  phase  of  the  disturbance  during  the 
THPS  experiment  from  the  results  calculated  by 
superposing  the  perturbations  generated  by  ilPS 
could  be  attributed  to  non-linear  effects.  An 
example  of  such  a  summation  is  given  in  figure 
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10  (for  F=104  at  R^*=1100  and  Z=0).  The 

d 

upper  part  of  tiiis  figure  presents  the  two  iso¬ 
lated  lii^  disturbances  participating  in  the  in¬ 
teraction  while  the  lower  part  shows  the  mea¬ 
sured  THPS  disturbance  in  comparison  to  the 
calculated  linear  superposition. 


Fig.  10  Two  experimental  UPS  eigen  functions, 
their  linear  superposition  and  the  THPS 
eigen  functions  for  the  same  case, 
F=104,  R^«b|I00  (Xj^SScm),  all  on 

Z=0. 
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Fig.  11  Comparison  between  THPS  and  the 
linear  superposition  of  the  two  HPS  at 
phase  shifts  of  A^O  and  )r/2. 


The  resemblance  between  the  reconstructed  and 
measured  amplitudes  and  phases  of  the  pertur¬ 
bations  is  satisfactory.  The  comparison  shown  in 
figure  11  is  applied  to  other  phase  shifts 
between  the  original  sources  with  equal  success. 
Some  of  the  observed  differences  may  stem 
from  the  imperfections  in  the  assumed  symme¬ 
try  of  the  laminar  boundary  layer. 


It  can  be  concluded,  on  the  basis  of  the  very 
good  agreement  observed  in  figure  10,  that  a 
linear  superposition  is  permissible  within  the 
range  of  parameters  considered.  The  validity  of 
this  technique  is  also  demonstrated  in  figure  II 
(for  F=>I04  and  at  R^8=990)  for  the  spanwise 

distribution  of  the  disturbance.  In  this  case,  a 
simple  reflection  technique,  with  respect  to  the 
plane  of  symmetry  (i.e  Z=0  in  fig.  2),  was  used. 
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Streaniwise  Interaction 


Having  established  the  validity  of  the  linear 
superposition,  between  two  experimentally  mea¬ 
sured  harmonic  wave  trains  introduced  at  the 

same  X  (i.e.  the  same  while  being  separ- 

6 

ated  along  the  span,  encouraged  us  to  expand 
the  numerical  scheme  to  include  other  possible 
spatial  interactions.  We  thus  considered  the  in¬ 
teraction  ol'  two  equal  frequency  (F=I04)  wave 
trains,  originating  from  point-sources  located  in 
tandem  (one  downstream  of  the  other)  at 
R  *=700  and  800.  It  should  be  emphasized  that 

O 

the  calculated  patterns  of  amplitude  and  phase 
arc  only  representing  the  Y  location  at  which 
the  amplitude  is  maximum  and  neglect  the 
changes  in  the  shape  of  the  eigen  functioits 
occurring  along  the  span.  The  validity  of  these 
results  is  also  limited  to  the  lower  amplitudes  of 
the  ensuing  disturbance,  thus  assuring  linear 
behavior.  The  calculated  amplitude  and  phase 
distributions  of  the  two  isolated  HPS  distur¬ 
bances  are  plotted  vs  Re2  in  Figure  12.  In  the 
chosen  example,  the  two  superimposed  wave 
trains  originated  from  point-sources  located  at 

R^*=700  and  800  respectively,  while  |he  result- 
0 

ing  disturbance  is  calculated  and  presented  at 

R^*=900. 

0 


Fig.  12  Calculated  amplitude  and  phase  distri¬ 
butions  for  two  isolated  MFS  originat¬ 
ing  at  Rf*=700  and  800  and  presented 

at  R^*=900. 


The  individual  amplitudes  of  each  source  were 
approximately  equal  on  the  centerline  (figure 
12)  at  Rr*=900  but  the  normalized  spanwise 

distribution  of  the  amplitude  depended  on  the 
distance  from  the  source  and  on  the  R^*  at  the 

source.  The  result  of  the  superposition  is  shown 
in  Figure  13  for  three  relative  phases  of  the  dis¬ 
turbance  on  the  centerline,  namely  A^O,  tr/2 
and  r.  It  is  clear  that  only  a  local  suppression  of 
the  harmonic  disturbance  is  feasible  in  spite  of 
the  fact  that  the  relative  choice  of  the  ampli¬ 
tudes  is  arbitrary.  When  complete  cancellation  is 
attained  on  the  plane  of  symmetry  the  result  at 
Re2  greater  than  2S,(X)0  is  not  affected  by  the 
interfering  wave.  One  may  continue  the  compu¬ 
tation  of  the  interacting  wave  trains  to  other 
streamwise  locations  (Fig.  14)  and  observe  that 
the  amplitude-peak  moves  away  from  the  plane 
of  symmetry  with  increasing  X. 


Fig.  13  Calculated  amplitudes  of  a  linear  sup¬ 
erposition  of  the  two  isolated  HPS  dis¬ 
turbances  presented  in  Fig.  12  at  three 
relative  centerline  phases. 


The  effectiveness  of  this  suppression  is  mea¬ 
sured  by  comparing  the  ensuing  amplitudes  to 

the  amplitudes  existing  in  the  flow  without  in¬ 
terference  (i.e.  amplitudes  produced  by  an 
isolated  harmonic  point-source  disturbance).  In 
the  present  example  an  attempt  was  made  to 
suppress  a  disturbance  originating  at  Rg*=700 

by  another  disturbance  located  directly  down¬ 
stream  at  Rj*=800,  while  the  results  were 

observed  at  various  spanwise  locations  between 
R^*=850-I050.  Even  the  most  effective  sup¬ 
pression  (occuring  at  A^x)  resulted  initially  (at 
R.«  <  l()00)  in  absolute  amplitudes  along  the 

Q 

span  which  were  higher  than  the  amplitude  pro¬ 
duced  by  a  single  HPS.  Further  downstream,  the 
interference  is  favorable  and  the  combined 
amplitude  of  the  perturbation  is  lower  than  for 
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Fig.  14  Calculated  amplitude  vs  Re2  of  stre- 

ainwise  THPS  interaction  (originating 

at  R^«»700  and  800)  plotted  at  various 
6 

downstream  (R=R^*)  stations. 
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Rez/1000 

Fig.  IS  Calculated  amplitude  vs  Re2  of  an 
isolated  HPS  interaction  (originating  at 
R^«=700)  plotted  at  various  down¬ 
stream  (RsR^e)  stations.  Symbols  sim¬ 
ilar  to  Fig.  14. 


an  isolated  HPS  (fig  IS).  The  rise  in  the  initial 
amplitude  of  the  combined  disturbance  at  short 
distances  from  the  secondary  source  may  intro¬ 
duce  non-linear  effects  which  were  not  consi¬ 
dered  presently  and  consequently  the  eventual 
suppression  predicted  for  R^*  >  1000  may  not 

materialize. 


Acknoledgements 


This  project  was  supported  in  part  by  a  grant 
from  the  United  States  Air 
Force  (AFOSR  89-0307)  and  monitored  by  Dr. 
J.  McMichael. 


Refrences 


Conclusions 

li  was  confirmed  experimentally  that  the  evolu¬ 
tion  of  a  single,  harmonic,  point-source  distur¬ 
bance,  in  a  boundary  layer  is  well  predicted  by 
a  linear  model.  Furthermore,  the  spatial  interac¬ 
tion  of  two  or  more  such  low  amplitude  distur¬ 
bances,  can  be  modelled  by  linear  superposition 
in  spite  of  the  complicated  phase  relationships 
existing  among  such  disturbances.  This  simplifi¬ 
cation  enables  one  to  predict  the  net  effect  of 
any  desired  distribution  of  point-source  distur¬ 
bances,  as  long  as  the  resulting  amplitudes 
remain  low. 

The  total  suppression  of  a  harmonic,  three 
dimensional,  disturbance  within  a  short  stream- 
wise  distance  from  its  detection  is  virtually 
impossible  to  achieve  unless  one  contemplates  a 
distributed  control  mechanism  consisting  of 
large  number  of  individually  activated  sources 
of  perturbation.  This  conclusion  should  also  be 
valid  for  the  impulsively  generated  three  dimen¬ 
sional  wave-packet  because  its  spanwise  struc¬ 
ture  consists  of  Fourier  components,  which  are 
individually  identical  to  the  harmonic-point- 
source  disturbance  (Ref.  14). 


1.  Kachanov,  YU.S.,  Development  of  spatial 
wave  packet  in  boundary  layer,  proceedings 
of  the  lUTAM  symposium  on  Laminar- 
turbulent  transition,  Novosibirsk  1984,  Ed. 
Kozlov,  V.V.,  p.  115. 

2.  Mack,  L.M.  and  Kendall,  J.M.,  Wave  pat¬ 
terns  produced  by  a  localized  harmonic 
source  in  a  Blasius  boundary  layer,  1983, 
AIAA  paper  83-(X)46. 

3.  Wehrmann,  US  patent  number  3,362,663, 
1968. 

4.  Liepmann,  H.W.,  Brown,  G.L.  and  Nosen- 
chack,  D.M.,  Control  of  laminar  -instability 
waves  using  a  new  technique,  J.  Fluid 
Mech.,  1982,  vol.  118,  p.  187. 

5.  Thomas,  A.S.W.,  The  control  of  boundary 
layer  transition  using  a  wave  superposition 
principle,  J.  Fluid  Mech.,  1983,  vol.  137,  p. 
233. 

6.  Pupator,  P.T.  and  Saric,  W.S.,  Control  of 
random  disturbances  in  a  Boundary  Layer, 
Abstract  submitted  to  the  2nd  shear  flow 
control  conference,  1989. 

7.  Caster,  M.,  A  Theoretical  Model  for  the 
development  of  a  wave  packet  in  a  laminar 

boundary- layers,  1975,  Proc.  Roy.  Soc. 
Lond.,  A.  347,  P.  271. 


7.12 


8.  Caster,  M.  and  Grant,  1.,  An  experimental 
investigation  of  the  formation  and  develop^ 
mem  of  a  wave  packet  in  a  laminar  boun¬ 
dary  layer,  1975.  Proc.  Roy.  Soc.  Lond.,  A. 
347,  P.  233. 

9.  Benjamin,  T.B.,  The  development  of  three 
dimensional  disturbances  in  an  unstable 
film  of  liquid  flowing  down  an  inclined 
plane,  J.  Fluid  Mech.,  10,  1961,  p.  401-419. 

10.  Criininale,  W.O.  and  Kovasnay,  L.S.G.,  he 
growth  of  localized  disturbances  in  a  lami¬ 
nar  boundary  layer,  J.  Fluid  Mech.,  1962, 
14,  p.  59. 

11.  Squire,  H.B.,  On  the  stability  of  three- 
dintensional  disturbances  of  viscous  flow 
between  parallel  walls,  Proc.  Roy.  Soc. 
Lond.,  1933,  A.  142,  p.  621. 

12.  Benny,  D.J.  and  Gustavsson,  L.H.,  A  new 
mechanism  for  linear  and  nonlinear  hydro- 
dynamic  instability.  Stud,  in  App.  Math., 
1981,  64:185-209. 

13.  Katz,  Y.,  On  the  evolution  of  a  turbulent 
spot  in  a  laminar  boundary  layer  with  a 
favorable  pressure  gradient,  PhD  thesis, 
1987,  Tel-Aviv  Uni. 

14.  Seifert,  A.,  On  the  interaction  gf  low  amp¬ 
litude  disturbances  emanating  from  discrete 
points  in  a  Blasius  boundary  layer,  PhD 
thesis,  1990,  Tel-Aviv  Uni. 


7.13 


APPENDIX  B 


Active  Control  of  Skin  Friction  and  Separation  in  Certain  Class  of 
Wall-Bounded  Flows 
by  I  V^gnanski 


ACTIVE  CONTROL  OF  SKIN  FRICTION  AND  SEPARATION 
IN  CERTAIN  CLASS  OF  WALL-BOUNDED  FLOWS 


presented  by: 
I.  Wygoanski 


Department  of  Aero^ace  Engineering  University  of  Arizona 

and 

Department  of  Fluid  Dynamics  &  Heat  Transfer 
Tel  Aviv  University 


ABSTRACT 

Experimental  investigations  of  large  coherent  structures  in  turbulent  shear  flows  bypassed 
the  boundary  layer  in  the  presence  of  a  strong,  adverse- pressure- gradient  and  the  wall- jet.  Both 
are  wall  bounded  flows  having  one  characteristic  in  common:  their  mean  velocity  profile  is 
inviscidly  unstable  to  two  dimensional  perturbations.  Thus  the  identification  of  the  large  coherent 
structures  with  the  predominant  instability  modes  might  be  extended  to  this  class  of  flows  and 
could  be  quantitatively  analyzed.  The  apparent  similarity  between  the  wall-jet  and  a  combination 
of  a  free  jet  and  a  boundary  layer  is  explored  was  and  the  relevance  of  the  solid  surface  to  the 
evolution  of  the  large  coherent  structures  was  assessed.  From  this  point  of  view,  the  boundary 
layer  in  a  strong  adverse  pressure  gradient,  is  regarded  as  a  wake  evolving  in  the  vicinity  of  a 
solid  surface.  In  both  flows  the  significance  of  the  outer  region  is  accentuated  while  the  no-  slip 
conditions  at  the  solid  surface  are  maintained. 

The  response  of  these  flows  to  external,  two  dimensional  excitation  is  currently  being 
investigated  in  an  attempt  to  further  the  understanding  of  the  interactions  between  the  inner  and 
the  outer  structures  in  a  turbulent  boundary  layer.  While  doing  so.  two  technologically  important 
effects  were  discovered.  The  mean  flow  in  the  wall  wake  remained  attached  in  spite  of  the  strong 
adverse  pressure  gradient  which  caused  separation  in  the  absence  of  the  excitation  while  the  skin 
friction  in  the  wall  jet  was  reduced  as  a  consequence  of  the  excitation. 

INTRODUCTION 

Investigations  of  large  coherent  structures  in  turbulent  shear  flows  concentrated  either  on 
wall  bounded  flows  like  a  boundary  layer  or  a  channel*’^  or  on  free  shear  flows  like  the  mixing 
layer,  the  wake  and  the  jet^'^.  In  the  latter  category  of  flows,  the  large  coherent  structures  were 
clearly  identified  as  the  predominant  insubility  modes  and  were  quantiutively  analyzed  in  this 
context.  In  fact,  one  may  formally  show  that  the  coherent  eddies  existing  in  a  two-dimensional, 
turbulent,  mean  flow  are,  to  a  first  approximation,  governed  by  the  Orr-Sommerfeld  equation 
provided  one  decomposes  the  equations  of  motion  into  coherent  and  incoherent  motion^,  one 
neglects  the  second  order  products  of  both  motions  and  one  restricts  the  divergence  of  the  mean 
flow.  The  uiple  decomposition,  therefore,  enables  one  to  link  the  hydrodynamic  stability  analysis 
with  quantitative  definitions  of  coherent  structures. 

However,  the  large  coherent  structures  in  wall  bounded  flows  are  much  more  complex  and 
more  difficult  to  identify  then  in  free  shear  flows.  Many  forms  of  such  structures  were  observed 
visually'^-*  but  no  consensus  was  reached  as  to  their  origin  and  their  precise  association  with  the 
enhancement  of  the  skin- friction  or  with  the  rate  of  growth  of  the  boundary  layer.  Since  Hussain 
and  Reynolds^' failed  to  excite  growing  two-dimensional  modes  in  a  turbulent  channel  flow,  it 
was  tenutively  concluded  that  plane  instability  modes  are  not  related  to  the  large  coherent 
structures  observed  in  channels  and  boundary  layers.  However,  instabilities  akin  to  the  ones 


responsible  for  the  generation  of  hair-pin  vortices  in  a  laminar  boundary  layer**,  may  still  be 
identified  with  the  creation  of  large  coherent  structures  in  fully  turbulent  wall  bounded  flows.  In 
fact  "sublayer  streaks"  were  recently  modelled  by  Landahl*^,  by  using  an  instability  concept 
associated  with  three  dimensional  intermittent  disturbances  growing  algebraically  in  the  direction 
of  streaming. 

Free  shear  flows  are  inviscidly  unstable,  while  boundary  layers  in  the  absence  of  an  adverse 
pressure  gradient  are  not  (figure  1).  There  is  enough  evidence*^  to  suggest  that  a  boundary  layer 
on  the  verge  of  separation  resembles  a  free  mixing  layer  and  it  responds  to  external  stimuli,  like 
a  mixing  layer.  The  receptivity  of  this  boundary  layer  is  attributed  to  an  inviscid  instability 
which  may  dwarf  any  other  form  of  instability  and  generate  large,  predominantly  two- 
dimension^,  coherent  structures.  The  wall  jet  is  inviscidly  unsuble  (figure  1)  in  its  outer  region 
and  may  thus  possess  large  coherent  structures  characteristic  of  a  plane  turbulent  jet.  The  wall  jet 
might  be  the  ideal  flow  configuration  for  resolving  the  intricate  interactions  between  the  outer 
and  the  inner  structures  in  a  turbulent  boundary  layer,  because  it  offers  a  larger  degree  of 
flexibility  and  controllability  of  flow  parameters  than  a  boundary  layer  does,  regardless  of 
pressure  gradient.  This  flexibility  stems  from  the  fact  that,  the  vorticity  in  the  outer  layer 
depends  on  the  added  momentum  flux  which  for  a  given  free  stream  and  jet  efflux  velocity 
depends  also  on  the  dimension  of  the  nozzle.  Thus  the  susceptibility  of  the  wall  jet  to  external 
perturbations  might  also  depend  on  the  ratio  between  the  wave-length  of  the  perturbation,  the 
width  of  the  nozzle  and  the  upsueam  boundary  layer  thickness.  Consequently  the  structure  of  the 
wall  jet  can  be  progressively  altered  and  the  importance  of  the  outer  vortical  layer  on  the  wall 
region,  can  be  evaluated. 

The  boundary  layer  in  a  strong  adverse  pressure  gradient,  may  be  regarded  as  a  wake 
adjacent  to  a  solid  surface  and  might  be  subjected  to  the  same  kind  of  stability  analysis.  In  both 
flows  the  significance  of  the  outer  region  is  accentuated  while  the  no- slip  conditions  at  the  solid 
surface  are  maintained.  The  vorticity  near  the  surface  of  the  wall-jet  has  an  opposite  direction  to 
the  vorticity  in  its  outer  layer  and  may  detract  from  its  strength.  On  the  other  hand,  the  vorticity 
near  the  surface  of  the  wall- wake  may  reinforce  the  vorticity  in  its  outer  region.  The 
superposition  of  an  external  stream  on  the  wall- jet  weakens  the  vorticity  in  the  outer  region  until 
it  resembles  a  turbulent  boundary  layer  in  the  limit. 

Although  one  has  less  flexibility  with  the  wall- wake  than  with  the  wall- jet  in  an  external 
stream,  the  understanding  of  the  former  flow  is  technologically  very  important.  In  particular,  the 
control  of  the  large  scale  eddies  in  a  turbulent  boundary  layer  having  a  vanishingly  small  wall 
stress*^  might  lead  to  a  delay  in  boundary  layer  separation  which  is  omnipresent  in  most  flows 
around  solid  bodies.  Preliminary  observations  associated  with  the  wall- jet  and  wall- wake  are 
presently  discussed,  stressing  some  features  common  to  both  flows.  Many  more  measurements 
might  be  required  before  the  concept  is  proven  but  the  results  obtained  to  date,  are  encouraging. 


THE  MEAN  FLOW  IN  THE  FORCED  WALL- WAKE 


The  flow  over  a  solid,  convex  wedge  was  investigated  (figure  2).  The  wedge  consisted  of 
two  long  flat  plates  which  were  hinged  togethti  to  provide  a  discontinuity  in  the  slope  of  the 
surface  and  thus  in  pressure  gradient.  An  18°  or  a  22°  wedge  angles,  which  could  not  normally 
have  been  negotiated  by  the  upstream  boundary  layer  without  separation,  were  used  in  the 
experiment.  The  different  wedge  angles  depended  on  the  detailed  geometrical  conditions  at  the 
discontinuity.  A  thin  metal  flap  spanned  the  entire  test  section  at  the  hinge  location  and.  in  the 
absence  of  the  forced  oscillations,  the  flow  separated  at  its  trailing  edge.  The  surface  of  the 
splitter  plate  was  roughend  upstream  in  order  to  ensure  that  the  boundary  layer  approaching  the 
discontinuity  was  is  thick,  turbulent  and  two-dimensional.  The  velocity  profiles  measured  are 
typical  for  a  turbulent  boundary  layer  in  this  range  of  Reynolds  numbers  (i.e.,  the  shape  factor  H 
a  5*/9=  1.6  for  Rej*  =  8(X)).  The  maximum  turbulent  intensity,  u'/U  is  approximately  7%  50 

mm  upstream  of  the  hinge.  The  displacement  amplitudes  of  the  flap  did  not  exceed  ±1  mm.  The 
imposed  streamwise  component  of  velocity  oscillations  were  measured  by  a  hot-wire  located  in 
the  potential  flow  near  the  flap.  The  velocity  perturbation  was  recorded  together  with  the  forcing 
signal,  and  the  data  were  ensemble  averaged  in  order  to  remove  the  free  stream  turbulence  of  the 
wind  tunnel  and  other  sources  of  random  disturbances  from  contaminating  the  results. 

The  flow  generated  downstream  of  separation  resembles  the  classical,  unforced  mixing  layer 
observed  in  the  absence  of  pressure  gradient.  The  vorticity  thickness,  obtained  by  fitting  a 
straight  line  to  the  mean  velocity  profile  at  the  point  of  inflection!^  increases  linearly  with  X,  as 
it  does  in  the  classical  mixing  layer.  The  slope  of  the  vorticity  thickness  dS^/dX  s  0.23  is  in 
good  agreement  with  previously  made  measurements  in  which  the  boundary  layer  on  the  splitter 
plate  was  tripped.^  The  pressure  recovery  coefficient  Cp  *  (P  -  Pjnh)/ j,  was  less  than  0.1 

over  the  distance  of  60  cm  because  the  flow  was  fully  separated  (figure  3). 

Activating  the  flap  at  frequencies  corresponding  to  Strouhal  numbers  based  on  initial 
momentum  thickness  (St  «  ^^hich  were  smaller  than  O.OI  and  reference  amplitudes 

>  0.25%,  resulted  in  a  reattachment  of  the  flow  to  the  surface  and  a  pressure  recovery  of 
Cp  =  0.5  over  the  corresponding  distance  considered  before  (figure  3).  The  flow  was  partially 
reattached  at  reference  amplitudes  as  low  as  0.1%.  One  may  conclude,  therefore,  that  the 
externally  introduced  oscillations  were  sufficiently  amplified  to  allow  the  boundary  layer  to 
overcome  the  adverse  pressure  gradient  imposed  by  the  divergence  of  the  walls  and  the  existence 
of  the  flap  at  the  discontinuity. 


Figure  3: 


The  p.e»i.e  diMribetion  "■  “f  "  C,  *i.h  « 

The  dependence  of  some  important  boundary  layer  parameters,  like  the  shape  factor  H,  and 
the  skin- friction  coefficient  Cj,  on  X  is  plotted  in  figure  4  for  some  of  the  frequencies  tested 
and  for  the  3  amplitudes  at  a  single  frequency  f  =  40  hr..  The  momentum  thickne.<!s  increases 
almost  linearly  with  X.  The  shape  factor  decreases  initially  with  X  and  then  gradually  increases 
towards  the  trailing  edge  of  the  plate,  while  the  skin  friction  coefficient  behaves  in  an  opposite 
fashion.  The  X  location  at  which  H  is  smallest  corresponds  to  the  location  at  which  Cf  is  largest, 
as  might  have  been  expected  from  the  momentum  integral  equation  provided  the  pressure 


increases  monoton ically  with  X.  The  location  at  which  H  has  a  minimum  appears  to  depend 
primarily  on  the  frequency  of  forcing  and  not  on  the  amplitude,  while  the  minimum  value  of  H 
attained  depends  mostly  on  the  amplitude  of  the  imposed  oscillations.  The  value  of  fX/U  at 
which  H  attains  its  minimum  is  approximately  2.  suggesting  that  the  wave-  length  of  the  imposed 
oscillations  determines  the  preferred  reattachment  distance. 


Figure  5: 


Mean  velocity  and  turbulent  intensity 
profiles  measured  at  x  =  10  cms 


U/Ut 


51 

u, 


(V.) 


Representative  mean  velocity  and  turbulent  intensity  profiles  measured  at  X  =  10  cm  are 
plotted  in  figure  S  for  three  frequencies  considered  at  a  forcing  amplitude  of  1%.  These  profiles 
are  typical  of  boundary  layers  at  moderately  adverse  pressure  gradients*^  with  the  exception  of 
the  maximum  level  of  turbulence  which  exceeds  20%  and  is  more  characteristic  of  the  level 
measured  in  separated  flows.'^  When  the  mean  velocity  profiles  are  plotted  in  the  usual  Taw-of- 
the-  wall"  coordinates  the  universal  logarithmic  distribution  remains  unchanged  regardless  of  the 
forcing  frequency  or  amplitude: 


U+  =  5.6  log  Y*  +  5  .  (1) 

The  maximum  value  of  Y'*'  at  which  the  equation  applies  depends  not  only  on  the  pressure 
gradient  but  also  on  the  amplitude  and  the  frequency  of  forcing.  It  can  be  best  correlated  with 
the  dependence  of  H  on  these  parameters.  The  general  representation  of  the  law- of- the- wake" 
also  applies  in  all  cases  considered,  i.e., 


W  =  C  sin^ 


(2) 


but  the  constant  C  depends  on  the  frequency  and  amplitude  of  the  forcing  (fig.  6).  It  seems  that 
in  a  "wall- wake"  there  is  a  large  region  in  which  the  turbulent  stress  is  approximately  constant  in 
contradistinction  to  the  measurements  done  in  a  wall  jet  at  moderate  Reynolds  numbers'^. 


Figure  6:  Mean  velocity  profiles  plotted  in  the  wall  coordinates. 


THE  TRANSFORMATION  OF  THE  WALL  JET  IN  STREAMING  FLOW 

TO  A  BOUNDARY  LAYER 


The  experimental  research  was  carried  out  in  a  low  turbulence  wind  tunnel.  The  wall-jet 
originated  from  a  two  dimensional  nozzle  flush  mounted  into  the  floor  of  the  tunnel  in  order  to 
minimize  its  effect  on  the  evolution  of  the  boundary  layer.  The  slot- width  was  altered  by  moving 
tiie  entire  front  plate  in  the  X- direction  and  inserting  appropriate  filling  strips  in  the  tunnel 
floor.  The  jet  turns  rapidly  in  the  direction  of  streaming  and  the  effects  of  its  deflection  could 
not  be  detected  beyond  10  slot  widths  in  most  test  conditions.  The  velocity  was  mostly  measured 
in  the  developed  region  of  the  wall  jet  starting  from  20  slot-  widths  downstream  of  the  nozzle 
and  extending  sometimes  beyond  200  slot- widths.  The  Reynolds  number  at  the  nozzle,  Rej,  was 
altered  by  changing  the  efflux  velocity,  Uj,  and  the  width  of  the  slot,  b.  The  ratio  between  the 
free  stream  velocity  and  the  jet  efflux  velocity  was  also  altered  by  varying  the  velocity  in  the 
tunnel,  ,  and  by  varying  Uj. 

Typical  stream  wise  evolution  of  mean  velocity  profiles  are  shown  in  fig.  7  where  the  top  set 
corresponds  to  an  initial  velocity  ratio  /Uj  =  0.085,  the  central  set  to  0.59  and  the  bottom  to 
0.93.  The  effect  of  the  upstream  boundary  layer  is  progressively  stronger  from  top  to  bottom. 
These  three  sets  of  data  were  chosen  because  they  may  represent  different  categories  of 
instability.  The  velocity  profiles  in  the  top  data-set  contain  only  one  inflection  point  in  the  outer 
region.  The  outer  region  in  this  case  resembles  a  free  jet  at  all  X  locations.  There  are  two 
inflection  points  in  the  outer  region  of  the  wail- jet  near  the  slot  of  the  second  set,  underlining 
the  importance  of  the  upstream  boundary  layer  in  this  case.  The  relatively  strong  Uj  overcomes 
this  velocity  defect  whereupon  downstream  of  the  fourth  station  the  mean- velocity  profile  no 
longer  contains  two  inflection  points.  The  effect  of  the  upstream  boundary  layer  is  much  stronger 
in  the  bottom  set  which  resembles  the  development  of  a  wake  in  the  vicinity  of  a  solid  surface. 
Since  the  ratio  of  the  free  stream  velocity  to  the  jet  velocity  approaches  unity  near  the  nozzle 
the  flow  appears  to  relax  directly  to  an  ordinary  boundary  layer  without  passing  through  the 
intermediate  stage  of  a  wall- jet  embedded  in  an  external  stream.  The  wall  jet  in  str  laming 
flow  may  therefore  represent  the  boandary  layer  in  the  limit. 
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Figure  7;  The  evolution  of  a  wail  jet  in  an  external  stream 


The  attainment  of  similar  mean  velocity  profiles  is  not  essential  in  analyzing  the  evolution  of 
the  large  coherent  structures  but  it  is  very  convenient  because  it  enables  one  to  reduce 
significantly  the  amount  of  computations  and  measurements  required  to  prove  the  concepts 
outlined.  It  is  well  known  that  a  complete  self  preservation  (i.e.  self  similarity  of  the  velocity 
disuibution  and  of  the  turbulent  stresses)  of  a  turbulent  wall-jet  may  only  be  possible  in 
especially  tailored  pressure  gradient‘s.  Nevertheless  many  investigators  managed  to  collapse  some 
of  the  data,  the  mean  velocity  profiles  in  particular,  on  self  similar  plots.  Using  a  Gallilean 
transformation  one  should  consider  (Umax-U^)  to  be  a  dominant  quantity  governing  the 
development  of  the  flow.  When  local  value  of  [(U-U.  )/(Umax-U„ )]  was  expressed  as  a 
function  fCq)  (where  q  ■  Y/Y„/2  represents  the  dimensionless  distance  from  the  wall  while 
is  the  distance  measured  from  the  surface  where  the  quantity  [(U-U. )/(Umax-U. )]  =  1/2) 
most  of  the  mean  velocity  profiles  measured  in  each  data-set,  collapsed  onto  a  single  curve 
(figure  8).  The  collapse  was  much  better  when  U«  /Uj  S  0.4  than  for  higher  ratios  of  U«  /Uj. 
The  shape  of  the  mean  velocity  profile  was  affected  by  the  ratio  U.  /Uj  while  a  change  in  Rej 
did  not  result  in  a  visible  alteration  of  the  normalized  shape  of  the  velocity  distribution.  The  rate 
of  spread  of  the  wall  jet  is  greatly  affected  by  the  velocity  ratio  and  by  Rej  (fig.9).  It  varies 
almost  linearly  in  the  direction  of  streaming  as  it  did  in  the  absence  of  external  flow. 


The  fully  developed  wall  jet  in  the  absence  of  an  external  stream  attains  a  local  equilibrium 
which  is  independent  of  the  detailed  conditions  at  the  nozzle‘^’^9  The  sole  parameter 
determining  the  evolution  of  the  jet  is  it's  initial,  kinematic  momentum  flux,  J.  This  approach 
scaled  correctly  the  three  most  important  parameters  in  the  wall  jet;  the  width  of  the  flow  Y^/2- 
the  decay  of  the  maximum  velocity  in  the  jet.  and  the  local  wall  stress,  x«. 

f 

:  ^-^-=F2(^)  :  Vp[-jJ  =F3(^)  (3) 

where: 

2 

^=[XJ/v^)  while  J*[Uj  b)  and  Uj  is  the  jet  velocity  at  the  nozzle  exit.  The  width  of  the  slot  is 
b;  p  and  v  are  the  density  and  kinematic  viscosity  of  the  fluid  respectively.  X  is  a  streamwise 
distance  measured  from  the  nozzle. 

It  made  these  parameters  independent  of  Rej  and  enabled  an  independent  evaluation  of 
The  question  of  scaling  took  different  proportions  than  originally  anticipated  because  it  appears 
that  the  flow  was  not  scaled  correctly  in  the  past.  The  addition  of  a  free  stream  complicates  the 
problem  because  one  regards  the  maximum  velocity  in  the  wall  jet  as  an  excess  velocity  over  the 


free  stream,  thus: 


Fl(0  =  [(U„  -  U«)vyj  (4) 

One  has  also  to  assume  that  the  paramount  parameter  governing  the  flow  is  the  excess 

r 

momentum  flux:  I  (u  -  U.  }udY  which  for  a  uniform  jet  velocity  in  the  absence  of  an 
•'o 

upstream  boundary  layer  degenerates  to  [b(Uj  -  U.  )Ujl  leading  to  two  possible  definitions  of 
the  streamwise  dimensionless  distance  4: 

4  ■  IXb(Uj  -  U«  )Uj  /v2]  or  ^  s  XIXb(Uj  -  U«  )Uj  /v^)  (5) 

where:  -  U«.  )/(Uj  +U_  ) 

When  the  upstream  boundary  layer  is  thick,  the  simplified  definition  of  excess  momentum  flux 
might  no  longer  be  adequate  because  one  has  to  account  for  the  momentum  loss  in  the  upstream 
boundary  layer. 

The  utility  of  the  definition  of  %  becomes  apparent  when  one  substitutes  the  observed,  self 
similar,  normalized  form  of  the  velocity  profile  (figure  8),  into  the  momentum  integral  equation 

which  makes  t^/p  (v/J)  a  function  of  X  and  on  [U.  /(Ug^-U^)]  while  being  independent  of 
RCj.  Therefore,  when  [U.  /(U„-U.)]  is  assumed  to  be  negligible,  the  similarity  conditions  of 
the  wall  jet  in  absence  of  an  external  stream  are  satisfied.  When  (Ug,  U. )  a  different  set  of 
conditions  for  self  similarity  prevail.  This  explains  some  additional  difficulties  in  analyzing  this 
flow. 


THE  FORCED  TURBULENT  WALL- JET  WITHOUT  AN  EXTERNAL  STREAM 

This  flow  was  examined  first  because  of  it  is  well  known  as  being  self  similar  with  respect  to  the 
selected  parameters.  External  exciution  of  free  shear  flows  had  a  strong  effect  on  a  variety  of 
time  averaged  quantities,  like  the  mean  velocity  distribution,  the  rate  of  spread  of  the  flow,  the 
turbulent  intensity  and  the  Reynolds  stress.  It  is  thus  natural  to  start  the  present  discussion  of  the 
experimental  results  by  exploring  parametrically  the  effects  of  forcing  on  the  same  quantities  in 
the  wall  jet. 

A  variety  of  forced  flows  were  considered  in  the  present  experiment,  they  differed  in  jet 
momentum;  slot  width;  frequency  and  amplitude  of  the  imposed  oscillations;  and  finally  the 
manner  in  which  forcing  was  introduced  to  the  flow.  A  change  in  each  one  of  the  above 
parameters  required  concomitant  traverses  across  the  flow  at  various  X  locations.  A  sample  of  the 
data  is  shown  in  figure  10.  The  velocity  profile  plotted  on  the  left  (figure  10a)  exhibits  self 
similarity  with  respect  to  Yg|/2  irrespective  of  the  X  location  while  the  profile  plotted  on 

the  right  (figure  1(%)  proves  the  existence  of  self  similarity  irrespective  of  the  changes  in  all 
other  parameters  considered.  Both  profiles  are  ctnnpared  to  the  dimensionless,  unforced,  velocity 
profile**  and  represented  by  a  solid  curves.  One  may  conclude  that  two  dimensional  external 
excitation  does  not  affect  the  normalized  form  of  the  velocity  distribution  within  the  range  of 
parameters  considered.  Replotting  the  velocities  adjacent  to  the  wall  on  an  expanded  scale  (figure 
lOc).  and  comparing  them  to  unforced  velocity  measurements  at  otherwise  identical  conditions, 
indicate  small  but  consistent  differences  between  the  two  sets  of  data.  These  differences  will  be 
discussed  later  in  some  detail. 
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Figure  10:  "  ^ 

Velocity  profiles  using  self-  similar  coordinates: 


We  were  surprised  to  discover  that  neither  the  rate  of  spread  nor  the  decay  of  the  velocity 
scale  were  affected  by  the  forcing.  The  solid  lines  drawn  in  figure  1 1  represent  the  unforced  data 
accumulated  for  Rej  2  7300  while  the  broken  lines  correspond  to  Rej  $  5000.  The  symbols 
plotted  represent  data  corresponding  to  various  frequencies  of  excitation  at  nozzle  Strouhal 
numbers  (Stj  ■  f*b/Uj)  varying  between  3  S  Stj*10*  S  17  (corresponding  to  St  ■  fv^/J^  of  0.4  S 
St*  10^^  i  444)  and  amplitudes,  based  on  the  stream  wise  velocity  perturbation  near  the  nozzle, 
ranging  from  2%  to  20%.  Although  some  consistent  differences  depending  on  the  level  of  forcing 
and  on  the  Strouhal  number  can  be  detected  they  do  not  alter  our  initial  conclusion  about  the 
lack  of  sensitivity  of  these  mean- flow  parameters  to  the  external  excitation.  The  cumulative 
effect  of  the  threshold  in  Rej  occurring  around  Rej  =  S(XX)  appears  to  be  much  more  significant 
then  that  of  the  two-dimensional  excitation.  Dimensional  analysis  of  the  independent  parameters 
in  the  externally  excited  wall- jet  suggests  that  variables  other  than  ^  might  affect  the  mean 
velocity  distribution.  Whenever  the  external  excitation  is  harmonic,  the  broadening  or  the 
distortion  of  mean  flow  will  depend  on  the  square  of  the  local  amplitude  of  the  velocity 
perturbations.  In  some  cases  the  finite  amplitude  might  be  introduced  directly  by  the  forcing 
mechanism,  while  in  others,  the  flow  itself  might  act  as  an  amplifier.  We  shall  now  examine  the 
global  effect  of  the  excitation  on  the  loss  of  momentum  resulting  from  skin- friction  drag.  The 
velocities  plotted  in  figures  10c  suggest  that  forcing  the  wall  jet  might  have  an  effect  on  the 
drag.  Since  the  local  skin  friction  was  most  conveniently  determined  from  the  slope  of  the  mean 
velocity  profile  near  the  wall**,  the  procedure  was  repeated  presently  in-  spite  of  the  fact  that  it 
is  not  considered  in  the  literature  as  being  very  reliable^* .  The  velocity  measurements  plotted  in 
figure  12  are  shown  in  dimensional  form  because  this  representation  enables  the  reader  to  assess 
the  quality  and  the  amount  of  data  acquired  in  the  viscous  sub- layer.  Sometimes  more  than  10 
data  points  were  included  in  the  linear  fit  made  to  the  velocity  profile  near  the  wall.  There  is  no 
doubt  that  even  a  relatively  low  level  of  forcing  (the  data  plotted  correspond  to  an  initial 
amplitude  of  2%)  results  in  a  reduction  of  the  wall  shear  stress,  t^.  Heat  loss  from  the  hot-wire 
to  the  wall  was  assumed  to  be  negligible^  since  the  flow  was  turbulent  and  most  of  the  data  was 
taken  at  distances  from  the  wall  ranging  from  SO  to  80  wire  diameters  while  the  Re  based  on  the 
wire  diameter  varied  from  0.7  to  4. 

Plotting  x^/p  (v/J)2  versus  ^  and  comparing  the  results  with  the  skin  friction  measured  in 
the  absence  of  forcing  indicates  clearly  the  differences  in  in  spite  of  the  logarithmic  scale 
chosen  (figure  13).  At  no  instance  did  %  increase  above  it's  nominal  unforced  value.  Reductions 
in  the  skin  friction  of  approximately  10%  were  prevalent  at  most  frequencies  corresponding  to 
initial  excitation  amplitudes  which  were  lower  or  equal  to  3%.  However,  reductions  in  of 
approximately  40%  were  also  recorded  by  forcing  at  much  higher  amplitudes  corresponding  to 


10%  or  20%  of  the  efflux  velocity  at  the  nozzle.  Alternately,  forcing  at  a  preselected  frequency 
which  is  amplified  by  the  flow  achieved  the  same  drag  reduction  at  a  much  lower  input 
amplitude  (figure  13). 

In  attempting  to  sort  out  the  independent  contribution  of  frequency;  amplitude  and  Reynolds 
number  on  the  local  one  may  assume  that  the  dimensionless  x^/p  (v/J)^  depends  on  the 
distance  from  the  nozzle  (XJ/v^)  and  on  the  local  amplitude  of  the  coherent  perturbations 
present  in  the  flow  regardless  of  their  origin.  We  computed  the  phase- locked,  ensemble- averaged 
velocity  signals  and  from  them  deduced  the  local  r.m.s.  levels  of  the  coherent  motion.  These 
quantities  were  then  integrated  across  the  flow  to  provide  a  measure  of  the  intensity  of  the 
coherent  motion  at  a  given  X  station.  It  transpired  that  a  local  coherent  motion  having  an  average 
amplitude  of  0.3%  may  be  responsible  for  a  reduction  in  of  approximately  15%  while  an 
additional  increase  in  local  amplitude  by  a  factor  of  10  resulted  in  an  incremental  reduction  in 
'^w/P  (v/J)^  of  an  additional  15%  only.  This  suggests  that  the  reduction  in  the  local  skin  friction 
is  not  linearly  dependent  on  the  strength  of  the  coherent  motion. 


A  comparison  was  made  between  the  mean  velocity  distributions  of  the  forced  and  the  unforced 
wall- jets  in  the  vicinity  of  the  surface  by  using  the  conventional  "wall  coordinates".  It  seems  that 
external  excitation  modifies  the  velocity  distribution  in  what  appears  to  be  the  logarithmic 
region.  It  was  shown  that  the  slope.  A,  of  the  logarithmic  profile  (U^  =  Alog  B)  appears  to 
be  a  universal  constant  and  is  equal  to  5.5  while  the  additive  constant  B  was  strongly  dependent 
on  Rej  (4  <  B  <  10).  This  conclusion  applies  to  the  forced  wall- jet  as  well.  It  also  appears  that 
the  logarithmic  fit  to  the  forced  data  may  only  be  applied  at  larger  values  of  Y'^  when  compared 
to  the  unexcited  flow. 


Figure  12: 

The  effect  of  forcing  on  dU/dY  near  the  wall. 
Figure  13: 

The  evolution  of  the  dimensionless  skin  friction 


with  ^ 


Some  of  the  underlying  assumptions  used  in  deriving  the  logarithmic  velocity  profile  do  not 
apply  to  the  wall  jet  (i.e.  that  neither  the  logarithmic  region  is  far  removed  from  the  location  at 
which  Y*Y„  nor  is  the  stress  constant  within  that  region).  Consequently  the  existence  of  the 
logarithmic  region  appears  to  be  somewhat  fortuitous.  The  introduction  of  forcing  did  not  only 
affect  %  but  also  reduced  the  extent  of  the  constant  stress  layer  (fig.  14).  The  total  and  the 
viscous  stress  distributions  in  the  inner  part  of  the  wail- jet  are  plotted  at  the  top  of  this  figure  in 
order  to  assess  the  effects  of  excitation  on  these  quantities.  The  wall- stress,  which  is  entirely 
viscous  at  the  surface,  is  reduced  by  the  two  dimensional  forcing.  However  for  Y*  >  10.  the 
effect  of  excitation  on  the  viscous  stress  vanishes  since  forcing  seizes  to  have  an  effect  on  the 


mean  velocity  profiles  in  this  region,  and  thus  on  v(dU/dY).  At  >  30  where  the  contribution 
of  the  viscous  stress  to  the  total  stress  becomes  vanishingly  small,  the  latter  decreases  rapidly  with 
increasing  distance  from  the  wall  and  changes  sign  between  60  <  Y^  <  100  depending  on  the 
forcing  and  on  Re. 

The  lateral  distribution  of  turbulent  energy  production  [-(u'v^dU/dY]  with  increasing  distance 
from  the  surface,  also  plotted  on  figure  14,  indicates  that  external  excitation  reduces  the  the 
Reynolds  stress  and  with  it  the  turbulent  production.  The  integrated  turbulent  production  in  the 
inner  region  of  the  wall- jet  between  the  surface  and  the  location  at  which  Y  =  Y„  defined  by: 


PRODUCTION*  (6) 


is  reduced  by  an  approximate  factor  of  3  when  the  flow  was  forced  at  St  10*^  =  S.S  at  X/b  =  70 
(i.e.  forcing  reduced  the  value  of  the  production  integral  from  0.668  to  0.214).  This  is  a 
significant  effect  which  presumably  alters  the  entire  turbulent  energy  balance  in  this  flow  and 
not  Just  the  scales  of  the  large  coherent  eddies.  Since  the  mean  velocity  gradient  vanishes  at  Y'*^  = 
ISO,  the  Y  locations  at  which  u'v'  =  0  and  at  which  dU/dY  *  0  do  not  coincide^  leading  to  a 
region  of  a  weak  but  negative  turbulent  production  just  below  the  location  at  which  Y=Y„. 


Forcing  the  wall-jet  reduces  the  turbulent  intensity  in  the  vicinity  of  the  surface  and  this 
reduction  is  strongly  affected  by  St  while  all  other  parameters  are  maintained  constant.  The 
reduction  in  u'/Uj  exceeds  in  some  places  25%  as  may  be  deduced  from  figure  IS  (see  arrow 
comparing  one  dau  point  at  X/b*40  and  St*S.5  which  provides  an  example).  Since  the 
production  of  turbulence  increases  the  level  of  u',  the  reduction  in  u'  due  to  forcing  might  be  a 


direct  result  of  the  reduced  production  level  near  the  wall. 


Figure  14: 

The  distribution 
of  total  and 
viscous  stress 
and  turbulent 
production  near 
the  surface. 


Figure  15: 

The  distribution  of  u'  near  the  surface 
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THE  EVOLUTION  OF  HARMONIC  PERTURBATIONS  IN  THE  WALL-JET  AND 
WALL-WAKE 


The  propagation  and  amplification  of  small  amplitude  wavy  disturbances,  in  slightly  divergent 
free  shear  flows  was  discussed  by  Crighton  and  Gaster^^  and  others.  The  application  of  this 
analysis  to  the  wall  jet  is,  in  principle  the  same,  with  the  exception  of  the  'Vio  slip"  boundary 
conditions  at  the  wall  and  the  need  to  use  the  viscous  form  of  the  Orr- Sommerfeld  equation. 


This  problem  too  was  analyzed  by  Caster^  in  conjunction  with  the  stability  of  the  divergent 
boundary  layer.  The  stability  of  the  non- divergent,  laminar  wall  jet  was  considered  by  Tsuji  et 
ai^.  We  shall  therefore,  not  repeat  these  derivations  and  only  dwell  on  some  specific  aspects 
stemming  from  the  application  of  this  model  to  the  turbulent  wall  jet. 

Harmonic  excitation  of  both  flows  enabled  us  to  digitize  the  measured  velocities  at  any 
streamwise  location  together  with  the  forcing  signal  used  to  deflect  the  flap  located  at  the  apex 
of  the  wedge  in  the  wall-  wake  case  or  near  the  nozzle  in  the  wall  jet.  This  signal,  plotted  at  the 
bottom  of  each  figure  (in  fig.  16)  provided  the  necessary  phase  reference  with  respect  to  which 
other  measurements  taken  anywhere  in  the  flow  could  be  referred  and  ensemble- averaged.  Since 
there  is  a  possibility  of  a  subharmonic  resonance  the  long  time  series  were  subdivided  into 
segments  which  contain  two  waves  of  the  forcing  frequency.  Samples  of  such  phase-  locked  and 
ensemble  averaged  traces  of  the  streamwise  component  of  velocity  are  plotted  above  the  forcing 
signal.  The  periodicity  of  the  velocity  signal  is  seif  evident  but  its  harmonic  distortion  varies  with 
Y.  The  least  distorted  signals  appear  at  the  outer  edge  of  the  flow  where  the  velocity  is  not 
contaminated  by  turbulence. 

One  may  observe  a  slow  phase  advance  in  the  maximum  velocity  recorded  in  the  wall  jet 
with  increasing  distance  from  the  wall  provided  Y/Y„^2  ^  ^  location,  which 

approximately  coincides  with  the  inflection  point  of  the  mean  velocity  profile,  a  sudden  phase 
shift  occurs  which  is  followed  by  a  gradual  phase  delay  with  increasing  Y.  The  data  presented 
for  the  wall- wake  accentuates  the  sudden  phase  shift  in  the  outer  flow.  One  may  Fourier 
decompose  these  signals  and  plot  the  transverse  distributions  of  amplitude  and  phase  of  the 
signals  at  the  frequency  of  forcing.  The  transverse  distribution  thus  obtained,  is  compared  with 
the  amplitude  distributions  calculated  from  the  linear  subility  model.  Both  measured  and 
calculated  amplitude  distributions  are  normalized  by  their  local  maxima  which  occurs  fairly  close 
to  the  wall. 
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Figure  16:  The  amplitude  distribution  in  a  forced  wall  jet  -comparison  with  linear 
model.  A  the  amplitude  and  phase  distribution  in  a  wall-  wake 


When  similar  comparison  was  done  in  highly  unstable,  unbounded  turbulent  shear  flows  like 
the  mixing  layer,  the  interaction  between  the  coherent  motion  and  the  incoherent  turbulence  was 
not  accounted  for.  Since  these  flows  were  dominated  by  an  inviscid  instability  the  inviscid  Orr- 
Somerfeld  equation  was  solved,  and  the  role  of  the  Reynolds  number  did  not  have  to  be 
considered^.  The  primary  instability  in  these  flows  is  affected  by  viscosity  in  qjite  of  the  fact 
that  they  are  also  inviscidly  unstable.  The  no- slip  condition  at  the  surface  provides  significant 
viscous  stresses  and  enhances  the  dissipation  which  may  cause  the  harmonic  motion  to  decay.  The 
incoherent  turbulent  fluctuations  may  increase  this  dissipation  and  the  non  linear  interactions 
between  the  coherent  and  incoherent  fluctuations  may  provide  a  cascade  mechanism  through 
which  the  effects  of  viscosity  are  enhanced.  All  these  interactions,  which  are  not  accounted  for 
in  the  model,  might  be  lumped  together  into  an  equivalent  viscous  term  by  introducing  an  eddy 
viscosity^- 29.  Since  the  Reynolds  number  appears  as  a  parameter  in  the  Orr- Somerfeld  equation 
a  choice  of  a  fictitious  Re  introduces  an  indeterminacy  not  present  in  the  inviscid  calculations. 
For  example,  the  experimental  results  plotted  in  figure  16  agree  reasonably  well  with  theoretical 


model  provided  the  eddy  viscosity  was  10  times  larger  than  the  viscosity  of  the  fluid. 

The  disturbances  at  the  fundamental  frequency  in  a  wall- wake  may  also  be  viewed  as  an 
array  of  spanwise  vortices  extending  across  the  entire  boundary  layer  (fig.l6).  The  ISO**  shift  in 
phase  which  depicts  the  center  of  the  eddies  is  convected  along  the  surface  at  approximately 
one- half  of  the  free  stream  velocity  while  moving  away  from  the  surface  with  increasing  X. 
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Figure  17:  Spanwise  coherence  in  a  wall-jet 
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We  shall  now  explore  the  changes  that  took  place  in  the  large  scale  structures  in  response  to 
the  external  excitation.  We  shall  particularly  look  for  effects  occurring  near  the  solid  surface  in 
order  to  ascertain  the  reasons  for  the  reduction  in  x^.  The  improvement  in  the  two- 
dimensionality  of  the  large  scale  structures  is  one  such  effect  which  had  been  observed  in  both 
the  mixing  layer  and  in  the  wake^-^.  The  coherence  spectra  of  the  velocity  fluctuations,  sensed 
by  two  probes  separated  in  the  spanwise  direction  provides  a  convenient  measure  of  the  degree 
of  two  dimensionality  attributable  to  each  scale  (or  frequency)  of  the  motion.  Two  point 
correlations  of  the  broad- frequency  signal  also  yield  a  measure  of  two- dimensionality  which  is 
biased  somewhat  towards  the  most  energetic  eddies.  The  data  plotted  in  figure  17  represents 
coherence  at  spanwise  separation  distances  ranging  from  10  to  130mm  corresponding  typically  to 
more  than  S  characteristic  widths  of  the  flow.  The  measurements  were  done  at  a  variety  of  Y 
locations  in  the  presence  and  in  the  absence  of  external  excitation.  The  frequency  chosen  to  be 
shown  for  the  wall  jet  was  34  Hz  corresponding  to  Stsfv3/J^=88*10''^.  The  results  plotted  in 
figure  17  correspond  to  four  Y  locations:  starting  at  the  outer  edge  of  the  wall- jet  (Y/Y^/j  = 
proceeding  to  the  vicinity  of  the  inflection  point  of  the  mean  velocity  profile  (Y/Y„/2  =  ^)> 
going  down  to  the  location  at  which  U  =  and  ending  in  the  vicinity  of  the  wall.  The 
spanwise  coherence  in  the  outer  part  of  the  unforced  jet  is  negligible  provided  AZ/Y„/2  (^'8- 

I7a).  For  ^ZIY^i  ^  2  the  coherence  in  the  unforced  flow  exceeds  the  value  of  O.IS  over  a 
broad  range  of  low  frequencies  and  Y  locations,  suggesting  that  the  large  scale  structures 
occurring  naturally  in  the  flow  are  somewhat  coherent  over  this  span.  There  is  also  a  distinct 
peak  in  the  coherence  level  around  30  Hz  within  the  inner  part  of  the  unforced  wall  jet  which  is 
even  noticeable  at  AZ/Y,„/2  *3.  If  this  St  is  associated  with  the  natural  coherence  of  the  large 
eddies,  then  the  forcing  frequency  used  in  the  experiment  is  very  close  to  being  the  frequency 
which  undergoes  naturally,  the  highest  level  of  amplification  at  this  particular  X  location. 

External  exciution  increased  the  coherence  to  0.7S  at  the  outer  edge  of  the  wall- jet  (fig. 


17a)  but  failed  to  have  any  effect  near  Y/Y„/2  “  *  (^*8*  i7b).  This  effect  may  be  explained  with 
the  aid  of  figure  16  which  is  representative  of  the  coherent  amplitude  distribution  of  the 
streamwise  velocity  fluctuations  in  the  forced  wall- jet.  The  maximum  amplitude  of  the  coherent 
motion  in  the  outer  part  of  the  jet  occurs  around  Y/Y„^2  while  the  minimum  occurs  at 
Y/Y„/2  »  1  and  therefore  one  should  not  expect  any  improvement  in  the  coherence  in  the 
vicinity  of  this  Y  location.  The  improvement  in  the  coherence  observed  in  the  inner  region  of  the 
wall- jet  (at  yi^tu/2  *  dt  0.01)  can  also  be  associated  with  the  relatively  large  amplitude  of 
the  local  oscillations.  It  is  interesting  to  note  that  the  q;)anwise  coherence  in  the  forced  wall- jet 
did  not  diminish  by  increasing  AZ  and  it  was  actually  higher  near  the  surface  than  in  the  outer 
edge  of  the  jet. 

Similar  measurements  are  also  done  in  the  wall  wake  (fig.18)  where  the  maximum  coherence  in 
the  unforced  flow  is  approximately  0.3  beyond  a  separation  distance  of  AZ  =  30mm  which  is 
commensurate  with  the  local  boundary  layer  thickness.  External  excitation  increased  this  value  to 
0.8  over  most  of  the  boundary  layer  with  the  exception  of  the  height  above  the  wall  where  the 
local  velocity  U/U|  s  0.7.  This  height  corresponds  to  the  average  location  of  the  centers  of  the 
passing  eddies  as  they  were  detected  by  the  180°  phase  shift  in  the  phase  locked  amplitudes. 
Three  dimensional  plots  of  the  ^anwise  correlation  coefficient  taken  at  various  heights  above  the 
surface  are  shown  in  figure  19  for  the  forced  and  unforced  flows.  The  contour  plots  representing 
the  natural  flow  exhibit  some  periodic  behavior  near  the  surface  of  the  wail  jet  but  none  near  the 
surface  of  the  wall  wake  because  of  the  massive  separation  which  occurs  under  these  conditions. 
Forcing,  not  only  increased  the  spanwise  extent  of  the  highly  correlated  region  but  made  it  much 
more  periodic  in  time.  The  primary  eddies  appear  to  be  two-dimensional  and  regular. 


•  n  tt  m  m  m  tm  m 
AZ  (mol  AZ  liMi) 

Figure  18:  Maximum  spanwise  coherence  in  a  wall- wake 
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Figure  19:  Correlation  measurements  in  the  T-Z  plane  of  a  wall- wake 

The  unexpected  reduction  in  the  mean- skin- friction  prompted  us  to  examine  the  phase  locked 
distribution  of  this  quantity  over  two  periods  of  the  forcing  frequency.  The  data  was  based  on 
the  phase  locked  analysis  of  [3(U)/3Y]^  measured  at  various  streamwise  locations.  The 
oscillations  in  the  skin  friction  corresponded  to  the  oscillations  imposed  on  the  flow.  It  is 
interesting  to  note  that  the  maximum  skin  friction  observed  during  any  phase  of  the  forced 
oscillation  at  any  streamwise  location  was  approximately  equal  to  the  mean  skin  friction  of  the 
unforced  flow. 


CONCLUSIONS 


Two  dimensional  excitation  of  the  wall- jet  in  the  absence  of  an  external  stream  has  no 
appreciable  effect  on  the  rate  of  ^read  of  the  jet  nw  on  the  decay  of  it's  maximum  velocity. 
Careful  examination  of  the  flow  near  the  surface  (i.e.  for  0  <  Y'*'  <  100)  reveals  some  profound 
differences  which  manifest  themselves  in  reducing  the  skin  friction.  Local  reductions  of  30%  in 
the  wall  stress,  as  a  consequence  of  such  an  excitation,  were  not  uncommon.  The  skin  friction 
drag  which  is  the  only  contributor  to  the  loss  of  momentum  in  this  flow  was  also  reduced  by  a 
comparable  amount.  The  production  of  turbulent  energy  near  the  surface  was  reduced  resulting 
in  lowering  the  intensities  of  the  streamwise  component  of  the  velocity  fluctuations.  These 
effects,  which  were  observed  in  the  fully  developed  region  of  the  wall- jet  (i.e.  at  X/b  >  30).  are 
insensitive  to  the  method  of  fencing  but  they  are  sensitive  to  the  frequency  and  the  amplitude  of 
the  excitation. 

External  excitation  enhanced  the  two- dimensionality  and  the  periodicity  of  the  coherent  motion. 
This  enhancement  is  clearly  visible  near  the  surface  and  near  the  free  interface  of  the  turbulent 
flow.  The  large  coherent  structures  in  this  flow  might  be  identified  with  the  most  amplified, 
primary  instability  modes,  of  the  mean  velocity  profile.  Detailed  stability  analysis  confirms  this 
proposition  though  not  at  the  same  level  of  accuracy  as  it  did  in  many  free-  shear  flows. 

The  preliminary  observations  in  the  wall  wake  lead  to  similar  conclusions  and  represent  also  a 
successful  attempt  to  conuol  and  delay  the  separation  of  a  turbulent  boundary  layer.  The 
introduction  of  harmonic,  two-dimensional  oscillations  results  in  a  reattachment  of  the  flow  and 
changes  the  proportions  between  the  "wake"  and  the  "wall"  functions  whose  linear  combination 
represents  the  streamwise  velocity  distribution  in  a  turbulent  boundary  layer.  It  does  not, 
however,  alter  the  universal  form  of  these  functiems. 
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